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The Larmor precession of a neutral spinning particle in a magnetic field 

<D . 

\ confined to the region of a one dimensional-rectangular barrier is investigated 

for both a nonrelativistic and a relativistic incoming particle. The spin pre- 
cession serves as a clock to measure the time spent by a quantum particle 
traversing a potential barrier. With the help of general spin coherent state it 
\ is explicitly shown that the precession time is equal to the dwell time in both 

the nonrelativistic and relativistic cases. We also present a numerical estima- 

Q_i| tion of the precession time showing an apparent superluminal tunneling. 
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I. INTRODUCTION 

The time spent by a quantum particle passing through a potential barrier has been one 
of the most controversial question since the founding of quantum physics and has attracted 
considerable attention from both a theoretical perspective and a experimental view [l| 0. 
There are various approaches dealing with tunneling time Q , but there has been no clear-cut 
answer to this old question [[| 0. Recently a number of experiments [[]] |7j indicating 
superluminal transmission of photons has renewed interest in this subject. 
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In quantum mechanics, time enters as a parameter rather than an observable. Thus, 
there is no direct way to calculate tunneling time. For particles with given energy there 
exist at least three different notions of the traversing time in the literature || [|Kj ||11|| , 
i.e., the Wigner-time, the Biittiker-Landauer time and the Larmor-time, corresponding to 
different criteria. The Wigner time accounts for how long it takes for the peak of a wave 
packet to emerge from the exit face of the tunnel barrier relative to the time that the peak 
of the incident wave packet arriving at the entrance face. The calculation of the Wigner 
time is based on an asymptotic treatment of tunneling as a scattering problem utilizing the 
method of stationary phase to calculate the position of the peak of a wave packet. This 
tunneling time is simply the derivative of the phase of the tunneling amplitude with respect 
to the energy of the particle. Buttiker and Landauer consider the case that the height of 
the barrier or the amplitude of the incident wave is modulated sinusoidally in time. They 
have found that if the frequency of the modulation is very low, the tunneling particle will 
see the instantaneous height of the barrier and the transmitted waves adiabatically follow 
the modulation. However, as the frequency of the modulation increases, the transmitted 
waves will no longer be able to follow adiabatically the rapidly varying modulation. The 
Biittiker-Landauer time is the modulation period such that the transmitted wave begins to 
depart from an adiabatic following of the modulation. 

Larmor precession was first introduced long ago as a thought experiment designed to 
measure the time associated with scattering events Subsequently the method was 



applied to measure the tunneling time of particles penetrating barrier with a magnetic field 
confined to the barrier region, causing the spin of particle to precess [|L3|1 . The original 



scheme jD^] considered only the rotation of the spin in the plane that is perpendicular to the 
magnetic field. Later it was recognized that a particle tunneling through a barrier in the 
magnetic field does not actually perform a simple Larmor precession in a plane JTT]. The 



main effect of the magnetic field is to align the spin with the field since the particle with 
spin parallel to the magnetic field has lower energy and less decay rate in barrier region than 
that of particle with spin antiparallel to the magnetic field. The total angular change of 



the tunneling particle divided by the Larmor precession frequency is the Larmor time ||TT| 



The literature also invokes a dwell time [HJ defined as the ratio of integrated probability 
density over the barrier region to the incident flux. The dwell time measures how long the 
matter wave in the barrier regardless of whether the particle is reflected or transmitted. 

In the present paper we revisit the Larmor precession of a neutral spinning particle with 
a general spin coherent state as a clock to measure the tunneling time through a barrier 
and extend for the first time the study of quantum tunneling to the relativistic regime. The 
advantage of using the spin coherent state is that an equation of motion for the expectation 
value of spin operator in the magnetic field within a barrier is obtained and identified with 
the equation of spin precession. With the help of the equation of motion we find that both 
nonrelativistic and relativistic neutral spin 1/2 particles perform a simple Larmor precession 
in three-dimensional space. The Lamor precession time of spin in a magnetic field confined 
to the potential barrier is compared with the time in the absence of barrier in order to show 
the apparent superluminal tunneling. 



II. TUNNELING TIME FOR A NONRELATIVISTIC PARTICLE 



Consider the one-dimensional rectangular potential penetration for a neutral particle of 
spin | with momentum p and mass m. The Hamiltonian is 
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where Vo is the height of the barrier situated between -d and d , lol = is the Larmor 
frequency and /i, h denote the magnetic moment and Planck constant, respectively. Here 
the time-independent magnetic field B is assumed in the z-direction and is confined to the 
barrier. a± (72,(73 are the Pauli spin matrices. The incoming wave is in the x-direction 
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where kq = y denotes the wave number and E is the energy of the particle. The 
component U\ (u 2 ) of the incoming wavefunction corresponds to spin up (down). We assume 
the incoming spinor is a normalized spin coherent state, which is the eigenstate of spin 
operator o ■ n with unit eigenvalue [|H)|. 
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where n =(sin 9 cosy?, sin 9 sin <p, cos 9) is an arbitrary unit vector with a polar angle 9 and 
azimuthal angle ip. The two components of the spinor are found to be 
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We only consider the case of E < Vq for quantum tunneling . The wave function to the left 
of the barrier (x < —d) is 



and the transmitted wave function to the right of the barrier (x > d) is 
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In the barrier the wave function is 
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where K\ k 2 are given by 
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The coefficients Aj , B; , Q and Dj (i = 1,2) in the wavefunction are obtained from 



boundary conditions ipi(-d) = ^ 2 {-d) , rp 2 {d) = ^{d) and £ipi(x)\ x= - d = £ip2(x)\ 



\x=— d 



= iKo + i^ eidKo - dKiDi 

2 k,; 
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where 



T 



Ri 



(/^ + k 2 ) 2 sinh 2 (2ck;) + 4k§k? 

(k 2 , + /t 2 ) 2 sinh 2 (2cki) 
^ 2 + K2) 2 sinh 2 (2^)+4K 2 K 2 



<j>i = arctan( — l - tanh(2d^)) (10) 

2KQKi 

For our purposes we shall consider the case of infinitesimal field limit, such that 



«i^«-2^. ^ + ^' -=^ {V °- E) (U) 

The transmission and reflection probabilities can be expanded as the power series of the 
small quantity The first order approximation is 

rp rp, X rp f \ 9T TTIUJ L 8T TTIUJ L 

R l = R {Kl) ^R iK) ---^, R2 = R{K2) ^ R(K) + __± ( 12) 
and we have 



V7^~T(/t), *Jr^R 2 ~R(k) (13) 

If we denote the transmission and reflection probabilities of spin up (down) by T + (TL) 
and i? + (-R-) respectively, then it can be easily shown that 
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which indicates the conservation of probability. From the viewpoint of scattering, the out- 
going wave packet consists of both a reflected and a transmitted wave packets, which are 
separated from each other. The outgoing wave packet must be normalized to unity, since 
the incoming wave packet is normalized to unity. 

The expectation values of spin for the transmitted wave ip t = (^) in the infinitesimal 
field limit are 

/c\ _ hT(. K ) sin#cos(0 2 - <t>i + V?) 
{b 1 ) t - 



(S2) t 



2 n K )-^^cos0 

h T(k) sin 9 sin(0 2 — <t>i + V 9 ) 
2 T(k)-^*£ cos* 



Wt 1T(.)-^4cos9 (16) 

The expectation values of spin for the reflected wave ip r = (^) in the infinitesimal field 
limit are 

, . _ HR(k) sin6'cos(02 - <j>i + V?) 

,,T " 5 K«)-¥S«»' 

/c \ _ ^-R(^) sin^sin(0 2 - 0i + v 9 ) 

The second terms in the expectation values of the z-component of spin for both reflected 
and transmitted waves have an obvious interpretation as given in Ref.[ll] that the particle 
with spin parallel to the magnetic field has lower energy and less decay rate in barrier region 
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than that of particle with spin antiparallel to the magnetic field. Equations (16) and (17) 
show that the spin still performs a Larmor precession around the z-axis which cannot be 
realized using the special spin-polarization of an incoming particle with polar angle 6 — | 
and azimuthal angle p = 0. To see the spin precession explicitly we may take the sum 
of expectation values of spin components for the reflected and transmitted waves with an 
infinitesimal magnetic field, i.e. ($) = ^^^^^ = (A |S| A) + {A |S| A)- 
We have 

(51) = ^ sin 9 cos(0 2 -<i>i + (p) 

h 

(5 2 ) = -sin6 sin(0 2 - (f>i + <p) 

(5 3 ) =^cos9 (18) 

which are of the same form as for the Larmor precession of a spin in a uniform magnetic 
field. To see this let us consider a neutral particle in a uniform constant magnetic field B 
along the z-direction in the absence of potential barrier. The Larmor precession is obtained 
by solving the Heisenberg equation 

jS(t) =±[8(t),H a ] (19) 



with the spin Hamiltonian 



H s = ~hw L a 3 (20) 



If the initial wave function is given by the spin coherent state i.e. 
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the expectation values of the spin components at time t are 

h 

(Si(t)) = — sin 6 cos(— uj^t + ip) 
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(S 2 (t)) = ^sin6sin(-uj L t + ip) 

(S 3 (t)) =^cos9 (22) 



Comparing Eqs.(18) and Eqs.(22) the Larmor tunneling time tl is obviously obtained as 

T L = LU^ifa - <j> 2 ) (23) 

Using Eq.(ll) we expand in Eq.(lO) up to the first order of the small quantity 
. the Larmor tunneling time is found to be 

ttikq 4(1k(k 2 — Kq) + (k 2 + Kq) sinh(4d/«)] 
Tl ~ ~h~K Ak 2 k 2 + (k 2 + k 2 ) 2 sinh 2 (2ck) ^ ' 

We assume that the incoming particle is a neutron with energy E and the width and 

o 

height of the rectangular barrier are 2d=8A , Vo=470MeV respectively. The Larmor tunnel- 
ing time as a function of the particle energy E is shown in Fig. 1(a). The peculiar feature, 
however, characteristic of Larmor tunneling time is that it increases with the energy of 



incoming particle monotonically in agreement with the observation in Ref. [IB]. It is inter- 
esting to compare the Larmor tunneling time with the Larmor time of a neutron traversing 
a constant magnetic field B confined in region —d < x < d , but without a barrier. With the 
same procedure as that for the case with a barrier we find that the transmission probability 
tends to one in the small field limit. The Larmor time of passage through the magnetic field 
region in the absence of a barrier is 

which is just the ratio of the traveling distance 2d of the spinning particle to its speed 
v = — . Using the same parameter as in Fig. 1(a), the plot of r° as a function of energy E 
is shown in Fig. 1(b). The ratio r = p- can be smaller than one. In other word the speed of 
a neutron in a barrier is larger than that in the free space |T7[. For the parameters chosen 
here, if the speed of incoming particle is one tenth of the speed of the light in the vacuum, 
the speed of particle though the barrier would be 5.7 x 10 15 m/s. 

The dwell time is defined as the ratio of the probability P& of finding a particle within 
the barrier to the incident probability flux Jj 

r d = ^ (26) 



The incident probability flux Jj is 

= (27) 
m 

and the probability for the particle to be in the barrier is 

p fc = / i^ti^mdx 
J-d 

A(Ik(k 2 — «q) + (k 2 + Kg) sinh(4dK)] 
~ "k" 4K§K 2 + (K 2 + K§) 2 sinh 2 (2dK) 

The dwell time is found to coincide with the Larmor time exactly, = tl in agreement with 



(28) 



the result of Ref. |TT[ for the spin polarization perpendicular to the magnetic field. The 
identity = ti has also been demonstrated in Ref. |18[ for potential barrier of arbitrary 
shape. 

III. QUANTUM TUNNELING OF RELATIVISTIC PARTICLE 

A relativistic neutral particle of spin ~ with mass m and magnetic moment //, moving in 
an external electromagnetic field denoted by the field strength tensor F^, is described by a 
four- component spinor wave function ip obeying the Dirac-Pauli equation 

eft 1 

[Y-d, + mc 2 + -{jut^F^ = (29) 
where c is the velocity of light in vacuum, 7 M = (7°,7) are Dirac matrices satisfying 

{y\7"} = V (30) 

with g^=diag(l, —1, —1, —1), and 
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It can be shown that 



[r",7l (3i) 



\^ U F^ U = ia ■ E - £ • B (32) 



where E and B are the external electric and magnetic fields, a = 7%, (3 = 7 . Here we 
make use of the Pauli representation 
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yO (X; y 



i = 1,2,3 



(33) 



and Si = |Ej is the spin operator. 

We again consider the one-dimensional rectangular potential barrier. The magnetic field 
B assumed in the z-direction is confined within the barrier region. The Hamiltonian is seen 
to be 



Hd = caip x + (3mc 2 \x\ > d 
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When E<V , the wave function satisfying the stationary Dirac-Pauli equation 
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where 



and 
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The incoming wave, i.e. the first term on the right hand side of Eq.(36), is assumed 
to be a normalized spin coherent state as in Eq.(4). The coefficients Aj, Bj, Cj and 
(i = 1, 2) in the wavefunction are obtained from boundary conditions ipi(-d) = ip 2 (—d) and 
ip 2 (d) = foid), namely, 



where 
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To our purpose we again consider the infinitesimal field limit 
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where 



k = -J(mc 2 + V ) 2 -E 2 , C=-(mc 2 + V + E) (44) 
c c 

is the zero order approximation. The transmission and reflection probabilities can be ex- 
panded as the power series of the small quantity E £%£ . The first order approximation is 

Ri = m)*m-§%£, R 2 = R(k 2) ^R(k ) + §^ (45) 

and we have 

^/T\f 2 ~T(k), ^/R^R 2 -R{k) (46) 

(l + fo)(T(k) + R(k)) = l (47) 
The transition and reflection probabilities of spin up and down are given respectively by 

t + = (i + f^D^ = (i + tim^i 2 

= (I + f )\D 2 f = (I + f )T 2 \u 2 \ 2 
R + = (l + fS)\A 1 \ i = (l + fS)R 1 \u 1 \ 2 

R. = (1 + / 2 )|A 2 | 2 = (1 + f 2 )R 2 \u 2 \ 2 (48) 
We have probability conservation so 

T+ + R+ = \ Ul \ 2 , T„ + R„ = \u 2 \ 2 (49) 

T + + R + + T_ + R_ = 1 (50) 

The expectation values of spin for the transmitted wave are obtained in the infinitesimal 
field limit as 
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The sum of expectation values of spin components for the reflected and transmitted waves 
with an infinitesimal magnetic field is 

h 



(Si 



sin 9 cos ( 



0i + v) 



(S 2 ) = ^^-Jl S in 9 sm(fa- fa + p) 
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which are formally the same as the Larmor precession equation of spin £ in a magnetic field. 
To see this we solve the Heisenberg equation (19) with the Hamiltonian 



and the initial wave function 
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The expectation values of the spin components at time t is 

(Si(t)) = ^ sin 6 cos(— u L t + ip) 



(5 S (0> =^^p^oosfl (56) 



Using the approximation Eq.(43) for ki in Eq.(39) and Eq.(42), the Larmor tunneling 
time which is defined by tl = <t>1 ~^ 2 is obtained as 



/o AdktjE(k 2 ~ ft?) + Kck 2 + E£)(k 2 + f^) smMjr) 



For a relativistic neutron and a rectangular potential barrier of width 2d=8A and height 
Vo=6000MeV, the Larmor tunneling time as a function of kinetic energy Eh, which in 
relativistic case is defined as the total energy E minus the static energy mc 2 , is shown in 
Fig. 2(a) which is similar to the nonrelativistic case (Fig. 1(a)) except the time scale. The 
Larmor time through the magnetic field region in the absence of a barrier is 

which also is exactly the ratio of the traveling distance 2d to speed v = c\J 1 — (i^) 2 . 
Using the same parameters as Fig. 2(a), r° the function of the particle kinetic energy E k 
is plotted in Fig. 2(b). For the parameters chosen here, if the speed of incoming particle is 
2.9 x 10 8 m/s, the speed of particle tunneling through the barrier would be 6.4 x 10 15 m/s. It 
is interesting to see the difference numerically between the non-relativistic expression (24) 
and the relativistic formula (61). To this end, we plot the results from both expressions in 
Fig. (3) with the same kinetic energy of the incoming particles. 
In the Dirac theory, the incident probability flux Jj is 

Ji = ^caiipi = 2 ^° 2 (59) 

1 + Jo 

and the probability for the particle to be in the barrier is 
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P 6 = / ^m^mdx 
J-d 

/2[-4dfc(fc 2 - e)(k 2 - ffc 2 ) + Kk 2 + e 2 )(A: 2 + fie) smMjg)] 
(1 + / 2 ) Wo^ 2 + (A: 2 + / 2 e 2 ) 2 sinh 2 (2* )] 

The dwell time r d is the ratio of Eq. (60) to Eq.(59), 

Td 2cA; 4/ 2 A; 2 e 2 + (k 2 + ffc 2 ) 2 sinh 2 (^f ) 

Using the relation 

cj 2 - cA; 2 

it is obvious that the dwell time Eq.(61) equals exactly the Larmor time Eq.(57) 

IV. CONCLUSION 

Using spin coherent state of an incoming particle we show that a neutral spinning particle 
penetrating the potential barrier with a constant magnetic field gives rise to a Larmor 
precession from which the barrier interaction time i.e. a time length for particle to remain 
in the barrier is determined. The Larmor time coincides with the dwell time in both the 
nonrelativistic and relativistic cases. The numerical calculation shows that the Larmor 
tunneling time can be much smaller than the time that the particle penetrates a constant 
magnetic field without a barrier, which implies the apparent superluminal tunneling. 
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Figure Caption: 

Fig. 1(a) Larmor time calculated from the non-relativistic expression Eq.(24) in a barrier as a 

o 

function of the particle energy with 2d = 8A, m =1.67x 10~ 27 kg and V = 470MeV. 

Fig. 1(b) Larmor time of Eq.(25) without a barrier as a function of the particle energy with 
2d=8A and m = 1.67x 10~ 27 kg. 

Fig.2(a) Larmor time for the relativistic case Eq.(57) in a barrier as a function of the particle 
kinetic energy with 2d = 8A, m = 1.67xl0~ 27 kg and V = 6000MeV. 

Fig.2(b) Larmor time of Eq.(58) without a barrier as a function of the particle kinetic energy 
with 2d=8A and m = 1.67xl0~ 27 kg. 
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Fig.3 Larmor times as a function of the kinetic energy for both the non-relativistic (dotted line) 

o 

and relativistic (solid line) particles with the same barrier width 2d=8A but various heights: 
(a) V =1000MeV, (b) V =3000MeV, (c) V =8000MeV. 
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